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ABSTRACT 

Post-launch boresight of the IBEX-Lo instrument onboard the Interstellar Boundary Explorer 
(IBEX) is determined based on IBEX-Lo Star Sensor observations. Accurate information on the 
boresight of the neutral gas camera is essential for precise determination of interstellar gas flow 
parameters. Utilizing spin-phase information from the spacecraft attitude control system (ACS), 
positions of stars observed by the Star Sensor during two years of IBEX measurements were 
analyzed and compared with positions obtained from a star catalog. No statistically significant 
differences were observed beyond those expected from the pre-launch uncertainty in the Star 
Sensor mounting. Based on the star observations and their positions in the spacecraft reference 
system, pointing of the IBEX satellite spin axis was determined and compared with the pointing 
obtained from the ACS. Again, no statistically significant deviations were observed. We con- 
clude that no systematic correction for boresight geometry is needed in the analysis of IBEX-Lo 
observations to determine neutral interstellar gas fiow properties. A stack- up of uncertainties 
in attitude knowledge shows that the instantaneous IBEX-Lo pointing is determined to within 
~ 0.1° in both spin angle and elevation using either the Star Sensor or the ACS. Further, the 
Star Sensor can be used to independently determine the spacecraft spin axis. Thus, Star Sensor 
data can be used rehably to correct the spin phase when the Star Tracker (used by the ACS) is 
disabled by bright objects in its field-of-view. The Star Sensor can also determine the spin axis 
during most orbits and thus provides redundancy for the Star Tracker. 
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1. Introduction 

One of the key scientific objectives of the 
Interstellar Boundary E xplorer (IBEX) mission 
(jMcComas et al.ll2QQ9bl ) is the accurate determi- 
nation of the interstellar neutral gas flow direction 
during the key observation period from December 
through March each year with the IBEX-Lo sensor 
([Mobius et alJ l2QQ9') in order to deduce the inter- 
stellar (IS) flow vector outside the heliosphere us- 
ing gravitational deflection of the flow by the Sun. 
To arrive at results that are at least commensu- 
rate in measurement accuracy with the previous 
determination of the helium IS flow vector by 



Wittd (|2QQ4f ) the pointing direction of the IBEX- 
Lo sensor boresight must be determined with an 
accuracy better than 0.2° in all directions. 

Placing such a requirement on the overall toler- 



ance budget between the spacecraft attitude con- 
trol and the mounting of the IBEX-Lo sensor, in- 
cluding the entire stack-up of mechanical toler- 
ances and thermal settling of spacecraft compo- 
nents, would have demanded significant resources 
that were not available. Therefore, a relatively 
simple Star Sen sor was included in the IBEX-Lo 
sensor package ( Fuselier et al.|[2QQ9[ ). In this way, 
the infiow direction of the interstellar neutrals can 
be directly tied to the positions of stars in the sky 
as an absolute astronomical reference frame. In 
order to achieve the required co-alignment solely 
through careful mechanical assembly, the Star 
Sensor was mounted on a shared baseplate that 
also holds the IBEX-Lo collimator, which deter- 
mines the direction of the incoming neutral atoms. 

As a simple Star Sensor concept, we use an 
optical system with a combination of split-V and 
pinh ole aperture, which is often used for Sun sen- 
sors ("Puchon fc Ver mande 1980'; 'Ki m et aL 2005 : 
[Wang & Chun 2006). Because the Star Sensor is 
co-aligned with IBEX-Lo, it continuously scans 
an approximately 8°-wide strip (as determined by 
the Star Sensor aperture) along a great circle in 
the sky that is at 90° ± 7° from the Sun. Ah 
stars within this strip can potentially be used for 
position determination for a single orbit, during 
which the IBEX spin axis orie ntation is kept con- 
stant (iMcC omas et al.l ( 2009b ^: see also discussion 
in the further part of the paper). As described 
below, such a scheme allows us to determine the 
positions of stars in this visible strip from the rela- 
tive timing of the star signals through the two legs 
of the split-V aperture relative to the start of each 
spin, which yields information on the spin phase 
and the elevation (relative to the spin plane of the 
satellite). This information can then be related 
to the stars' position in an astronomical reference 
frame. Relating the pointing of the Star Sensor 
to astronomical coordinates aids in the analysis of 
data from the IBEX-Lo sensor and from the IBEX 
mission as a whole. 

First and foremost, comparing the observations 
of several stars during a few sample orbits estab- 
lishes the absolute pointing direction of the Star 
Sensor and of the IBEX-Lo boresight in the IBEX 
coordinate system. Therefore, the pointing of the 
sensor is known in absolute terms for all time pe- 
riods when the IBEX Attitude C ontrol System 
(ACS), including its Star Tracker (jScherrer et al 



2009[ ). is operational according to its specifica- 
tions. During intermittent outages of the ACS due 
to bhnding of the Star Tracker, which occur dur- 
ing periods when either Earth or the Moon get 
close to the Star Tracker field of view, the space- 
craft loses its synchronization with the astronom- 
ical coordinates and the spin phase starts to drift 
slowly. Because the mounting of the Star Sensor 
and Star Tracker do not change substantially in 
fiight, the Star Sensor signal can be used to "de- 
spin" the IBEX data accurately so that the start 
of each spin is exactly re-aligned with its nomi- 
nal occurrence, 3° before the IBEX-Lo boresight 
passes the southern ecliptic pole. Finally, for or- 
bits that bring at least two well identifiable stars 
into the Star Sensor field of view, the orientation 
of the spin axis can be determined independently 
of the IBEX ACS. Therefore, the IBEX-Lo Star 
Sensor provides partial redundancy for the IBEX 
ACS system. 

In this paper, we demonstrate how the ori- 
entation of the IBEX-Lo boresight in the IBEX 
coordinate system was determined, including its 
uncertainties relative to the Star Sensor and the 
IBEX satellite system. We show that the IBEX- 
Lo boresight direction is known during all IS fiow 
observations to a very high accuracy. We also 
demonstrate that a relatively simple Star Sensor 
can serve as a stand-alone attitude determination 
component for simple spacecraft missions. 

In section 2, we provide a functional descrip- 
tion of the Star Sensor, followed by the Star Sensor 
calibration results for the accuracy of the determi- 
nation of the pointing direction and its brightness 
sensitivity for a range of stellar magnitudes in Sec- 
tion 3. Section 4 contains a description of the Star 
Sensor data and of the methods used to analyze 
them. Section 5 discusses the accuracy of spin 
axis determination from ACS data and presents 
analysis of the Star Sensor observations, starting 
from the determination of observed star positions 
in the sky in relation to their expected positions 
based on the IBEX spacecraft attitude determina- 
tion. Then the spin axis orientation is determined 
independently from the observation of at least two 
separate stars during individual IBEX orbits and 
compared with the IBEX attitude determination. 
The paper ends with 2 Appendices, in which the 
Star Sensor Simulation Program is presented and 
deviations of the spin axis derived from Star Sen- 
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sor observations in a few individual orbits from the 
ACS-determined spin axis are plotted. 

2. Functional Description of the IBEX-Lo 
Star Sensor 

The Star Sensor is to provide true pointing 
accuracy for the neutral atom observations rela- 
tive to known star positions in the sky. This is 
achieved through the use of the timing of the ar- 
rival of starlight through a succession of a split- 
V and pinhole aperture. The general scheme for 
the operation of the Star Sensor is shown in Fig. 
[TJ Along with IBEX-Lo, the Star Sensor points 
radially at 90° relative to the spin axis. As the 
satellite spins, the light from a certain star that 
is within the narrow angular band swept out by 
the Star Sensor field of view (FoV) arrive at the 
detector only when the star, one of the two legs of 
the split-V, and the pinhole aperture with a pho- 
tomultiplier (PMT) behind it are aligned. During 
each of these alignment events with a succession 
of partial and full overlap of the two apertures, as 
indicated in the view from the PMT through the 
pinhole aperture in the center panel of Fig. [T] for 
3 positions, the illumination of the detector be- 
hind the pinhole aperture increases and decreases 
sharply. The resulting shape of the light curve 
is close to a triangle, as indicated in the bottom 
panel of Fig. [T] with the 3 positions from above 
inserted. 

The mean timing of the two consecutive signals 
relative to the start of each spacecraft spin pro- 
vides azimuth information in the spacecraft coor- 
dinate system, and the timing difference between 
the two successive signals in the short sequence 
of a double-pulse, created by star position in the 
split-V aperture, provides the elevation. In this 
way consecutive sky circles will overlap due to the 
re-orientation of the spin axis towards the Sun at 
perigee of each orbit. During the first 2.5 years of 
operation, the orbital period of IBEX was kept at 
7-8 days, which requires a re-orientation maneu- 
ver of typically 6.9°-7.9° in each orbit. The Star 
Sensor is mounted in the shadow of the IBEX-Lo 
collimator to shield the sensor from sunlight. 

As can be seen from Fig. [T]the time difference 
between the occurrences of the two maxima de- 
pends on the actual elevation of the star relative to 
the Star Sensor boresight. The arithmetic mean of 
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Fig. 1. — Schematic view of the Star Sensor oper- 
ation. While the satellite spins around the vertical 
axis, the light from a star reaches the PMT when 
either one of the legs of the split V is aligned with 
the pinhole aperture and the star, thus leading to 
two consecutive near triangular pulses. Top panel: 
Side view with a star aligned through the right leg. 
Center panel: View seen from the PMT through 
pinhole and split-V aperture for three positions of 
the split V aperture (1, 2, 3) during a spin. Bot- 
tom panel: Light curve elapsing from right to left 
with positions 1, 2, and 3 from the center panel 
representations marked. 



the maxima occurrences determines the crossing of 
the boresight of the Star Sensor in azimuth or spin 
angle. The light curve for each star consists of two 
triangular shaped peaks with a full-width at half 
maximum (FWHM) that is equivalent to the an- 
gular width of the aperture system. Fig. [2] shows 
the aperture layouts for the Star Sensor that has 
been tailored to its requirements. With an aper- 
ture diameter dA = 1 mm, a respective width of 
1 mm for the two legs of the split V, its height 
of 6 mm, and an optical tube length Is = 40 mm 
between the two apertures. The FWHM of each 
of the two signals is 1.45° and the FoV in eleva- 
tion from —5.0° to +3.50°. The open diameter 
to house the aperture system is 11 mm. The pin- 
hole aperture diameter of 1 mm translates into 
an aperture area of 0.0079 cm^, which transmits 
1.23 X 10~^^ W into the detector for a Sun-like star 
of magnitude 3. The detector system consists of a 
compact photomultiplier (PMT) and an amplifier 
chain provides for light sensitivity down to about 
magnitude 4. The amplifier chain has a time con- 
stant of 12 ms, which shifts the Star Sensor signal 
by 0.3° in spin phase, given a nominal spin rate of 
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Fig. 2. — Star Sensor scheme, aperture configura- 
tions, and relevant linear dimensions. 



4.2 rpm. 

The IBEX spin period is typically 14.3 ± 0.2 s, 
which is equivalent to ~ 4.2 rpm. The PMT signal 
is accumulated into an array of 720 angular bins 
with an integration time of ~ 20 ms for 0.5° incre- 
ments. While the array has 720 bins, their width 
is not exactly 0.5°, but generally a little wider, 
with the last bin of each spin at reduced width to 
complete 360°. The reason is that the bin width is 
determined by an integer number k of basic time 
ticks with a granularity of 69.4 /is, where k is de- 
termined by the actual IBEX spin period as ob- 
tained from the IBEX ACS. k is always rounded 
up to the next integer that makes the bin width 
just > 0.5° so that always < 720 bins are filled 
with data for one spin. The Star Sensor signal 
is then accumulated over 64 spins. The absolute 
reference for this array of 720 bins is the space- 
craft spin pulse, whose timing is accurate to 10 jas 
(and thus negligible for required angle uncertain- 
ties < 0.1°) and which defines the beginning edge 
of spin bin with high angular accuracy. The spin 
pulse, which is intended to connect the spacecraft 
spin phase with an astronomical reference system, 
is issued when the IBEX-Hi boresight passes —3° 
from the northern ecliptic pole (NEP). Conversely, 
the IBEX-Lo boresight and thus the Star Sensor 
boresight passes —3° from the southern ecliptic 
pole at this time. 

The Star Sensor data are transmitted as a 720- 
bin array of 8 bit values every 64 spins, or ap- 
proximately every 15 minutes with the IBEX-Lo 
housekeeping data. For the first six months of op- 



eration, these data represented actual accumula- 
tion over 64 spins. Upon the realization that the 
IBEX ACS would lose its reference each time when 
a bright object (Earth or Moon) enters the FoV 
of the IBEX Star Tracker and the orientation of 
the spin pulse would drift slowly, a change in the 
Star Sensor data accumulation was implemented. 
Because the Star Sensor signals are generally so 
intense that PMT noise is almost negligible, no 
accumulation of the signal over consecutive spins 
is necessary. Therefore, after the first six months 
of operation (upon completion of the first set of 



IBEX EN A maps (|McComas et al.l l2009al )). the 



data are taken as 1-spin samples for every 64-spin 
accumulation period. 

3. Calibration of the IBEX-Lo Star Sensor 

Prior to integration into the IBEX-Lo sensor, 
the Star Sensor was fully tested and calibrated. 
Two functionally identical models were built and 
calibrated, the fiight model and the fiight spare 
model. Both models showed almost identical char- 
acteristics, and, in particular, their pointing, as 
is relevant to this investigation, turned out to 
be identical and repeatable within the established 
pointing accuracy when compared for identical 
settings during calibrations. Here, we report the 
calibration results for the fiight model (FM). Rel- 
evant to this investigation, the FM was calibrated 
for its response to light equivalent to set star mag- 
nitudes, for its angular response and boresight 
pointing in spin angle and elevation. 

3.1. Calibration setup for absolute bore- 
sight reference 

The Star Sensor calibrations were performed 
in a low light calibration facility of the magne- 
tospheric and ionospheric research laboratory of 
the University of New Hampshire, which contains 
an Integrating Sphere light source (by Sphere Op- 
tics) that is calibrated to NIST standards within a 
clean dark room. The optics is set up on an opti- 
cal bench that is isolated from the fioor. In order 
to use the Integrating Sphere, which produces ho- 
mogenous light over a 10 cm diameter opening, a 
cover plate with a 4 mm diameter pinhole centered 
on the opening was added to simulate individual 
stars. 

The boresight of the Star Sensor was designed 
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Fig. 3. — Schematic Star Sensor calibration setup. 
Bottom: Optical setup with a laser level in- 
serted between Integrating Sphere and Star Sensor 
mount, which contains a precision mirror for align- 
ment in the x-y plane. Center: Setup for align- 
ment in the x-z plane, with the mirror in the cen- 
ter height position (0°). Top: Optical setup with 
Star Sensor in the center height position. 

and fabricated as perpendicular to the baseplate of 
the IBEX-Lo sensor, which also defines the bore- 
sight of the IBEX-Lo collimator. Through care- 
ful choices of machining techniques the mounting 
of the Star Sensor within IBEX-Lo was controlled 
such that the collimator boresight is perpendicular 
to the baseplate to ±0.066°. 

In order to calibrate the pointing of the Star 
Sensor relative to an absolute boresight refer- 
ence that is perpendicular to its mounting plate 
a precision jig for the calibration facility was de- 
signed, whose perpendicularity could be adjusted 
and tested on the optical bench with a laser level. 
A schematic view of the test setup is shown in 
Fig. [3l The Star Sensor is mounted on a base- 
plate in one of three height positions (1, 2, 3) that 
are separated by Ah = 14 cm and with position 
2 such that the Star Sensor boresight points ex- 
actly at the artificial star aperture when correct 
alignment of the mount is achieved (shown in the 
top part of Fig. [3j). The separation between the 
artificial star aperture and the pinhole aperture of 
the Star Sensor is / = 248 cm, thus leading to el- 
evation pointing of +3.23 ± 0.02° (1), 0° (2), and 
— 3.23±0.02° (3). The jig is attached to a precision 
rotation table whose rotation axis represents the 
z-axis and is centered on the pinhole aperture of 
the Star Sensor so that the Star Sensor can be ro- 
tated relative to the artificial star to simulate the 
spacecraft spin. For exact alignment, the rotation 
table is set to 0°. 



To achieve alignment a laser level that emits 
laser light along the same axis to both sides is 
mounted between the Integrating Sphere and the 
Star Sensor mount (shown in the center and bot- 
tom parts of Fig. [3j). The laser beam adjusts itself 
to exactly level. The laser level is adjusted in y and 
z directions so that it is exactly centered on the 
artificial star aperture. A precision polished metal 
plate with 4 markings at equal distance from the 
center of the metal plate is placed into position 2 
of the Star Sensor mount. Now the mount is ad- 
justed in the y and z directions so that the laser 
beam is centered between the 4 center markings. 
Finally, the mount is adjusted in angle in the x-y 
and in the x-z plane so that the laser beam is ex- 
actly refiected back into the beam exit of the laser 
level. These alignments are achieved to ±0.5 mm 
so that an overall alignment accuracy and thus 
Star Sensor reference in perpendicularity to the 
mounting plate of < 0.018° was achieved. 

For consistency we used a mirror that was 
placed at the position of the polished plate which 
refiects the laser beam back on itself. Hence we 
double the travel distance of the laser beam, and 
thus, we doubled the accuracy. 

3.2. Calibration of angular response and 
boresight pointing 

To calibrate the Star Sensor for its angular 
response the orientation of the Star Sensor is 
changed in increments of 0.5° with the rotation 
table when the Star Sensor is mounted in either 
of the three height positions to allow for the sim- 
ulation of three different elevations. Fig. |4] shows 
the result from such a scan for the FM and the FS 
in comparison with an illumination equivalent to 
a magnitude 1.5 star. The measurements were ob- 
tained with engineering model electronics, which 
allows an amplifier output range from to 10 V. 
0° on the angle scale refers to exact alignment of 
the boresight with the line-of-sight to the artifi- 
cial star. The results in Fig. |4] demonstrate that 
the alignment of the Star Sensor was highly repro- 
ducible, even between different Star Sensor mod- 
els. The reproducibility was also achieved before 
and after vibration testing of the Star Sensor. 

For exact boresight pointing, the center po- 
sition in spin angle was verified to be accurate 
within ±0.065°. The separation of the legs of 
the split V aperture in the center position (2) 
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Fig. 4. — Angular scan with the Star Sensor 
(FHght Model in blue, Flight Spare in red) in ele- 
vation position 2 in 0.5° increments. 

in the mount, i.e. 0° elevation, was found as 
8.365 ± 0.039°. Including the fact that the two 
legs are inclined at 14.4° relative to the center line, 
the elevation pointing was found to be accurate to 
±0.077°. 

Including the mechanical tolerance stack-up of 
the collimator boresight of 0.066° and the uncer- 
tainty in the knowledge of the alignment of the 
Star Sensor in the calibration mount of 0.018°, 
the accuracies of the alignment of the Star Sensor 
and the collimator boresight are known to within 
0.094° in spin angle and 0.102° in elevation. 

3.3. Sensitivity and dynamic range of the 
Star Sensor 

Both Star Sensor models have been calibrated 
for a wide range of star magnitudes that include 
the magnitudes of the visible outer planets (Mars, 
Jupiter, and Saturn). This calibration provides 
the necessary information to set the bias voltage 
of the photomultiplier (PMT) tube adequately in 
flight and to help identify bright stars in the Star 
Sensor signals in flight. It should be noted that the 
star magnitudes are given as magnitudes in the vi- 
sual spectral band based on the color temperature 
of the lamp in the Integrating Sphere. To relate 
the results to actual star magnitudes the spectral 



curves of individual stars and the spectral curve of 
the PMT need to be factored in. Because the Star 
Sensor will be solely used for pointing purposes, 
no absolute calibration taking into account these 
spectral sensitivities has been attempted. The re- 
sult shows the wide dynamic range that the Star 
Sensor can operate in. In additional tests with 
a different light source, which was less well cali- 
brated than the Integrating Sphere, it was verified 
that the Star Sensor is even capable of reliable ob- 
servations of the half moon (magnitude ~ —10"^) 
with PMT voltages between 200 and 250 V, which 
are now used in flight for moon viewing time pe- 
riods. 

4. Data analysis 

4.1. Description of Star Sensor data 

The "first light" observations from the Star 
Sensor demonstrated that the instrument is work- 
ing excellently and very close to expectations. 
This is illustrated in Fig. [6] where an observa- 
tion from a selected orbit is compared to the sig- 
nal predicted by the Star Sensor Simulation Pro- 
gram (Appendix A). Several effects can obscure 
Star Sensor data: the signal from the Star Sen- 
sor can be saturated by the light of the Moon or 
Earth, when they are close to the FoV; and the 
synchronization of the spin pulse can be lost dur- 
ing a portion of the orbit if either the Earth or 
the Moon blind the spacecraft Star Tracker. The 
Star Sensor signals are used to correct pointing 
information when the spin-pulse is lost. 

Examples of such effects in the Star Sensor data 
are shown in Fig. [5l The upper two panels and 
the lower left panel show the Star Sensor signal 
as function of spin phase (vertical) and time (hor- 
izontal) in a color-coded representation, and the 
lower-right panel presents an example of the data 
from orbit 33 in a more detailed 3D view. 

In a scenario that is straightforward to analyze, 
a point source should show up as a pair of horizon- 
tal "rails", which correspond to the double peak 
of the star over time. Diffuse background should 
show up as broad horizontal bands. Any depar- 
tures from strictly horizontal structures mean that 
either a time-related phenomenon (like an object 
moving through the FoV) is observed or a problem 
with spin synchronization exists. The synchro- 
nization problems are illustrated in the upper left 
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Fig. 5. — Examples of data from various orbits, shown to illustrate the complex nature of the Star Sensor 
signal. The two upper and the lower left panels show a 2D time versus North Ecliptic Pole angle repre- 
sentation of the Star Sensor signal for 3 different orbits. The upper left panel shows an example when the 
star tracker is in eclipse, which causes the spin-phase to meander from its proper value. The lower-right 
panel presents a 3D visualization of the Star Sensor telemetry data collected from the entire orbit 33. The 
intensity of the signal is colour-coded and displayed at the vertical axis. See text for discussion. 



panel, where the horizontal structures are trans- 
formed into oblique wavy ones - and the waviness 
is similar for all elements of the picture. Only at 
the very end of this orbit, at the far right in the 
panel, the synchronization was restored. In such 
a case, the observations of point sources from the 
Star Sensor can potentially be used to despin the 
observations provided that at least one clearly vis- 
ible star is in the FoV. 

If the Star Sensor is to be used also to indepen- 
dently determine the spacecraft spin axis orienta- 
tion, at least two well separable stars are required. 
The upper right panel shows an example of a good 
orbit for this purpose: two stars are visible against 
a relatively weak background in the upper por- 
tion of the panel. There is even a possibility to 
resolve a third star, but this one is seen against 



a bright background strongly varying with spin 
phase, which makes using it for the spin axis point- 
ing determination problematic. However, the re- 
maining two stars are sufficient for the task. Only 
at the end of the orbit, strong reflected stray light 
from the Earth enters into detector and swamps 
the two good stars. 

The lower left panel illustrates a case of heavy 
contamination by the glow of the Earth (the big 
red structure in the upper portion of the panel) 
and Moon (the oval structure). There are also two 
bright stars shining through the background and 
the Moon glow. At about 2/3 into the orbit inter- 
val the Star Sensor was temporarily switched into 
the Moon Mode (with appropriate reduction of the 
PMT voltage) because of the expected passage of 
the Moon through the FoV. This time period is 
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visible as the vertical black strip. And indeed, 
the Moon signal was registered, as manifested by 
the two Hnes visible against the dark background 
cutting through the oval. The Moon shows a com- 
ponent of motion through the FoV from lower to 
higher elevation in the S/C reference system. As 
a consequence, the angular distance between the 
two peaks is increasing with time. During this 
particular orbit also the Earth was within FoV 
simultaneously with the Moon, as evidenced by 
the strong signal visible during the Moon Mode 
time period. The Earth signal is saturated even 
in the Moon Mode, nevertheless the Earth is re- 
solved and visible as a compact object (since it 
gives two distinct peaks). Because of the orbital 
motion of IBEX, the Earth is crossing the FoV 
downwards (i.e. its elevation decreasing, in con- 
trast to the Moon), as can be established from the 
fact the distance between the two peaks is decreas- 
ing with time. It is evident that the Earth and the 
Moon move also in azimuth, as indicated by the 
fact that the imaginary center lines between their 
double peaks are not horizontal, as are the sig- 
nals from the two stars. Since the synchronization 
is maintained in this orbit, the Star Sensor pro- 
vides us with time-resolved positional information 
on the objects moving through the FoV. Even un- 
der such challenging conditions, the signal from 
the two good stars at the beginning and the end 
of the orbit could be used for the axis pointing 
determination. 

This qualitative discussion is offered to justify 
that given the complex and varying composition of 
the data we decided not to develop an algorithm 
for an automatic selection of suitable data inter- 
vals for the analysis. Instead, we rely on human 
judgment supported by simulation. 

4.2. Determination of star positions in the 
spacecraft reference frame 

Determination of orientation of the optical axis 
of the IBEX-Lo instrument from the Star Sensor 
observations requires two steps: extraction of po- 
sitions of point-like objects (stars, planets) from 
the telemetry data and identification of these ob- 
jects with stars or planets listed in a stellar cat- 
alog or planetary ephemerides for the time of the 
observation. Typical reference objects most fre- 
quently observed by the Star Sensor will evidently 
be bright stars, but also the naked-eye outer plan- 



ets and even the Moon can be considered as addi- 
tional reference sources. 

The point-like stars are always observed against 
a comparably strong diffuse sky background. The 
most important components of the background are 
the zodiacal light, diffuse Galactic li ^ht and a huge 



number of faint, unresolved stars (jLeinert et al 



1998[ ). Hence the signal registered by the Star 



Sensor is composed of signals from the objects of 
interest for the Star Sensor on top of an extended 
background in the spectral sensitivity range of the 
detector. 

The method used to determine the positions of 
stars from the IBEX telemetry data is inverse to 
the method used to simulate the Star Sensor re- 
sponse function, described in Appendix A. The 
Star Sensor data collected during one entire orbit 
typically include a few hundred histograms (scan 
data blocks). To identify stars, we look for double 
pulses which are considerably stronger than the 
neighboring sky background level (Fig. [6]) and 
do not have other strong stars nearby. Strong 
neighbor stars in the FoV modify the shape of the 
pulses of the target star. Because the distortions 
of each of the two peaks usually differ, the dis- 
tance between the two peaks is changed, which 
impacts both the azimuth and elevation obtained 
from such a star. The selection of the candidate 
peak pairs is finalized through human involvement 
because of the varying and complex character of 
the signal, as discussed in the preceding section. 

Processing of observations from a given orbit 
begins with selection of a subset of data blocks 
free from synchronization issues and saturated re- 
fiected light. Then, the histogram portions with 
the peak pairs are determined. The latter choice 
is maintained for all suitable data blocks from the 
given orbit because the stability of the signal is 
very good for orbits in full spin synchronization. 
Only peaks stronger than a preselected level of 
0.15 V are analyzed. 

With the appropriate pairs of peaks selected, 
the local sky background level (approximated by 
straight line. Fig. [7|) is subtracted from the sensor 
signal, as illustrated in Figs [6] and [71 Based on the 
residual signal, the center position of each peak is 
determined through the calculation of its center- 
of-mass. Because the FWHM of the slit width is 
1.4°, each peak typically extends over seven an- 
gular bins. After identifying the maximum bin of 
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Fig. 6. — Illustration of the process of identification of the strongest vaHd peak pairs in the Star Sensor 
histogram. The Star Sensor telemetry histogram obtained during orbit 38 on 2009 Jul 20, 23:20:56 is shown 
as a red line. The simulated signal for the same orbit as obtained from the Star Sensor Simulator program is 
shown as a black line. The azimuth and brightness of stars brighter than m = 6 are shown (blue line). The 
azimuth of a star is the mean value of the position of two peaks that the Star Sensor returns from a point 
source in the sky. There are two bright stars with a high signal-to-background ratio in the figure. The star 
SAO 100944 is observed against a background whose intensity varies within the azimuth range occupied by 
the double peak. Thus, before processing of the star signal the background must be subtracted, which is 
done after fitting the background linearly (green line). 



a peak, the center-of-mass is obtained from the 
seven bins centered on this bin, using bin centers 
as angular coordinates. 

We obtain the spin angles ai, ol^ for both peaks 
of a pair. Before further processing of the selected 
pair we check whether it meets the criteria of the 
split-V aperture geometry, i.e. 



5.73° <OL2-OLi< 9.87° 



(1) 



If valid, the azimuth (spin angle) of the identified 
star astar in the spacecraft system is derived as 



Q^star = (<^1 + 0^2) /2 



(2) 



We do not require here that the peaks are equal 
in height because they still might be affected by 
unaccounted background and photocathode inho- 
mogeneities. 

With the azimuths of the two peaks and the 
azimuth of the star, the elevation of the star b^iar 
can be calculated from the formula: 



^star = arcsin 



tan ((^2 — OL\— Gy) /2) 



tan (/3y) 



(3) 



where = 14.4° is the tilt angle of each slit in 
the split-V aperture relative to the symmetry line 



and Gy = 8.4° is the angular separation of the 
slits at 0° elevation. 

4.3. Correction of an instrumental shift in 
spin angle (azimuth) 

A detailed comparison of the spin angles (az- 
imuths) of stars determined from the telemetry 
with the spin angles expected from simulation 
showed that the spin angles of the stars in the 
S/C frame are shifted relative to the positions ex- 
pected from simulations by an angle that increases 
with the spin angle. We identified two reasons for 
this effect: (1) a constant shift by ~ 0.3° is in- 
troduced by the Star Sensor amplifier, and (2) a 
linear "scaling" effect in the on-board signal for- 
mation process exists. The scaling effect comes 
up because of a finite resolution of the coding of 
the duration of the clock ticks used for bin width 
determination, indicated in Section 2. Here, we 
develop a correction factor, /g, that is the ratio of 
the actual bin width divided by the 0.5° bin width. 
The actual bin width can be solved for using a 
housekeeping register k that is used to solve for 
the time interval (in seconds) Ts = (/c + 192)/14400 
between the interrupts that separate each of the 
720 spin-bins. The actual bin-width (in degrees) is 
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Fig. 7. — Illustration of the method used for de- 
termination of the scanned star position in the 
S/C coordinates system from the Star Sensor his- 
togram data. In the telemetry signal (red line) 
there are two pulses, Al and A2, of the Star Sen- 
sor output voltage resulting from a star passage 
through the Star Sensor FoV. They are observed 
against sky background. With the boundaries of 
the peak pair in the histogram identified, the peaks 
are temporarily removed from the histogram, and 
a linear model is fitted to the remaining back- 
ground signal. Then the background is subtracted 
from the original signal, and the resulting peaks 
are processed to determine their positions and thus 
the star azimuth. The results are checked against 
the simulations. The simulation result is superim- 
posed on the data (marked as "simulation"). The 
position of the star along with all stars brighter 
than m = 6 is shown in the figure in blue. 

then A(j)s = SQO^Ts/Tcev^ where Tceu is the spin 
period of the S/C stored in the memory of the on- 
board computer and measured on a spin by spin 
basis. The correction factor fg = A<ps/0.5° or, 
combining factors: 



fs = {k + 192) 



720 



1 



14400 Tceu 



(4) 



The corrections for the instrumental shift (0.3°) 
and scahng of the spin phases are implemented 
directly in the program processing the Star Sensor 
telemetry data. 

The method for the star position determina- 
tion has been extensively tested on the signal ob- 
tained from simulations. A series of 50 cycles of 
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Fig. 8. — Illustration of the effectiveness of phase 
correction defined in Eq. HI Differences between 
the orbit-averaged observed and simulated az- 
imuth of observed before the correction are shown 
as the red histogram. The differences after the 
correction are represented by the blue histogram. 



the Star Sensor signal simulation/scanned star po- 
sition finding were performed for a few IBEX or- 
bits (spin axis positions). Since the simulations in- 
clude statistical fiuctuations of the intensity of the 
signal, the results were statistically distributed, 
but the mean position of the identified stars were 
found to be only slightly different from the actual 
star positions used for the simulation. For exam- 
ple, for IBEX Orbit 22 the error in right ascen- 
sion was only ~ 0.02° and in declination ~ 0.05°. 
In this way we made sure that the implemented 
method could determine the star positions with 
sufficient accuracy after the input data had been 
appropriately corrected. 

An illustration of quality of the correction de- 
fined in Eq. (j4]) is shown in Fig. [51 where his- 
tograms of the differences between the orbit- 
averaged observed and simulated spin angles (az- 
imuths) of the stars in the S/C frame are shown. 
The red histogram of pre-correction differences is 
not symmetrical and features an extended wing to- 
wards the negative differences. The differences be- 
tween the simulated and post-correction azimuths 
(blue in the figure) are much more symmetric 
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and the histogram is narrower. This histogram is 
shown in a finer scale in the right-hand panel of 

Fig. m 

4.4. Identification of found objects with 
catalog stars 

Once the positions of the stars in the spacecraft 
reference system are determined, they have to be 
identified with reference stars from the star cata- 
logue. The identification is based on a comparison 
of the telemetry data with simulations results (Fig. 
[6|) . The simulation of the signal expected from the 
telemetry is calculated for a given UTC time and 
spin axis position obtained from the Star Tracker 



signal processed by the ISOC (jSchwadron et al 
20091) . Based on the simulation results, a short 



list of candidate stars is returned whose positions 
are then compared with the observed star posi- 
tions. Details of the simulation program are given 
in Appendix A. Even though a software module 
to identify the stars from the observations was de- 
veloped and successfully tested, the manual iden- 
tification turned out to be more practical in sit- 
uations when the signal is affected by one of the 
unexpected complexities discussed in the preced- 
ing sections. 

4.5. Determination of the spin axis orien- 
tation 

If at least two stars for one orbit can be identi- 
fied with star catalog entries, it is possible to cal- 
culate the spin axis orientation. Determining the 
axis pointing requires finding a matrix R^* that 
transforms the orientation of the spacecraft refer- 
ence system, in which we have obtained the star 
positions, into the target inertial celestial refer- 
ence system. In our case, the target reference 
system is the equatorial system a nd the method 
adop ted is the TRIAD algorithm (jShuster fc Oh 
198ll ). There is relatively little robustness in the 
process of determining the S/C attitude using ob- 
servations of just 2 stars because there is no redun- 
dancy of information. Nonetheless, this method 
offers a very valuable way to compare the Star 
Sensor's attitude determination with that of the 
ACS system. Furthermore, the resulting spin axis 
pointings can be statistically analyzed to validate 
the actual pointing of the Star Sensor with that 
of IBEX-Lo in spacecraft coordinates. Finally, it 



is demonstrated that the Star Sensor could pro- 
vide the spacecraft attitude determination if such 
a need would arise. 

In order to determine the spin axis orientation 
with the use of the TRIAD algorithm, the posi- 
tions of the two stars must be determined quite 
accurately, which is satisfied by the selection cri- 
teria for the stars as described above. In practice, 
such conditions are not always fulfilled and thus 
there are orbits for which no determination of the 
spin axis position is possible because of the lack of 
suitable stars. 

The relation between the measured star posi- 
tion vectors in the S/C system vi^^ V2b and their 
position vectors in the equatorial system vu, V2i 
is: 

i;i5 = R^'vu and V2h = ^^'v2i (5) 

which must be solved for the transformation ma- 
trix R^* . The matrix algorithm is based on the con- 
struction of two triads and ti of ortho-normal 
unit vectors using the observed and reference star 
posit ion vectors Vu^ V2i^ 'Wib, V2b (Shuster fc OhI 
198lh : 



tb = [tib,t2b,tsi)] and ti = [tii,t2i,tsi] (6) 
where 



t2i 
tsi 

tlb 
t2b 
hb 

[^li, ^2i, ^3i] 
[^16,^26,^36] 



Vli 
Vli X V2i 
\Vii X V2i\ 

tli X t2i 

Vlb 
Vlb X V2b 
\Vib X V2b\ 

tlb X t2b 

= R'H = 

p6t| ^ 



(7) 



(8) 



pit 

R6t 



Finally: 



R 



bi 



= [^16, ^26, ^36] [^li, ^2i, ^3i]^ 



(9) 



(10) 



represents the transformation from the equatorial 
system into the S/C system. The inverse transfor- 
mation is the transposed matrix: 



R 



ib 



(11) 



11 



The IBEX spin axis is the z-axis of the S/C 
reference system. Therefore, the transformation 
matrix ahows the calculation of the spin axis Vzb = 
[0,0, 1] in the equatorial system as: 

t^« = (R''0^t^.6 (12) 

The method was implemented numerically and 
successfully applied to determine the spin axis 
pointing of IBEX during the first two years of sci- 
ence operations. 

5. Results 

5.1. Accuracy of spin axis pointing from 
ACS 

Pointing of the IBEX spin axis is determined 
by the ISOC based on data from the IBEX ACS. 
Experience from more than 2 years of operations 
showed that when the observation conditions are 
favorable, the ACS performs with high precision 
and according to specifications, enabling a very 
accurate determination of spin axis pointing. The 
spin axis is very stable during an entire orbit, with 
a very small precession amplitude of ~ 0.02°, i.e. 
comparable to the typical precision of the ACS 
system. The high stability of the IBEX spin axis 
enables us to adopt one solution for the entire du- 
ration of each orbit and to perform a reliable sta- 
tistical comparison with solutions obtained by the 
Star Sensor, as discussed in the following sections. 

5.2. Comparison of the predicted and 
measured positions of stars 

As a first verification step for the IBEX-Lo 
boresight in the spacecraft system, we analyzed 
the positions of all identified stars in the space- 
craft reference system for systematic deviations 
from the expected locations. For most of the suit- 
able orbits, positions of two stars could be deter- 
mined and the IBEX spin axis calculated. For a 
few orbits, the positions for either just one or as 
many as three stars could be determined. In total, 
69 stars were processed. The positions of stars in 
the S/C reference system were retrieved separately 
from each suitable telemetry block. Examples of 
the determinations of these positions are shown in 
Fig. M 

The positions of identified stars are somewhat 
distributed in both coordinates, with standard de- 



viations varying from orbit to orbit. An exam- 
ple of statistical fiuctuations of the differences be- 
tween the simulated and actually measured az- 
imuths and elevations for the star SA0246574 ob- 
served in Orbit 73 is shown in Fig. [TOl The po- 
sitions were stable and consistent. The median 
value of standard deviations in the phase angle 
was 0.02° and in elevation 0.06°. Such a behavior 
of a greater uncertainty in elevation than in spin 
angle was expected, as discussed in the previous 
section. 

The angular separation between stars found for 
a given orbit is very stable. The distribution of 
angular separations for an example for orbit 73 
has a standard deviation of less than ~ 0.04°. 

We calculated the mean positions of all identi- 
fied stars, followed by the mean deviations of these 
positions from the simulated values based on the 
spin axis positions obtained fro m the IBEX Sci- 
ence O perations Center (ISOC) Schwadron et al. 
(|2009[ ) using the spacecraft ACS. Histograms of 
the differences between the simulated and ob- 
served positions of the stars in the S/C reference 
system are shown in Fig. [TTl The uncertainty 
from the mean is given by ct/a/TV, where A" = 69 
is the number of observed stars. Therefore, the 
mean deviation in the phase angle is equal to 
-0.0155° ± 0.017°, while the mean deviation in 
elevation is -0.0179° ± 0.037°. About 68% of the 
results in elevation and about 85% results in spin 
angle are within ±0.2° from 0. We conclude that 
there is no statistically significant deviation of the 
Star Sensor boresight from the direction perpen- 
dicular to the ACS-determined spin axis to within 
0.037°, nor in spin angle relative to the spin pulse 
to within 0.017°. We conclude that there is no sta- 
tistically significant deviation of the FoV from the 
Star Sensor and thus the IBEX-Lo pointing other 
than expected from the observed defiection of the 
spin angle from spin axis and from the ground- 
measured uncertainty in the IBEX-Lo mounting. 
In addition, the Star Sensor signal can be used to 
correct the spin phase information with high pre- 
cision when bright objects, such as the Earth and 
Moon, bhnd the Star Tracker causing inaccuracies 
in the ACS data. 
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Fig. 9. — Examples of star positions in the spacecraft coordinate system for selected orbits. Black dots mark 
a position obtained from a single histogram block, the red error bars indicate the mean values with standard 
deviations. 



5.3. Independent spin axis pointing deter- 
mination using the Star Sensor 

For all orbits with at least two clearly identi- 
fied stars, we have calculated the pointing of the 
IBEX spin axis. We followed two paths: (1) we 
calculated the spin orientation separately for all 
spin blocks available from a given orbit, and (2) 
we calculated the spin axis from the mean star po- 
sitions averaged over a given orbit. Subsequently, 
we calculated the deviation of our pointing from 
the ACS pointing. The representation of the re- 
sults of the spin axis determination from Star Sen- 
sor observations are presented in the S/C reference 
system in Fig. [18] in Appendix B. The S/C refer- 
ence frame is oriented with the ACS-derived spin 
axis at the pole, and IBEX-Lo directed in the - 
Y direction. Typically, determinations for a given 
orbit form elongated streaks in the S/C reference 
system because the Star Sensor's accuracy for a 
star's elevation is lower than that of its azimuth. 
The quality of the IBEX pointing determinations 



varies between orbits because of the varying mea- 
surement conditions between orbits: relative posi- 
tion of the stars (differences in the azimuth and in 
the elevation), signal to background ratio of the 
stars, intensity and inhomogeneity of the back- 
ground, presence of the stray light, and the finite 
accuracy of the projection of the stars azimuth on 
the Star Sensor 720-bin histogram. The positions 
of the spin axis as obtained from individual data 
blocks for all processed orbits are shown in the left 
panel of Fig. [12] and the average positions for these 
orbits are shown in the right panel. The orbit av- 
eraged spin axis pointings determined by the Star 
Sensor are almost uniformly distributed around 
the pole. The small void region between 60° and 
330° azimuth is due to the season during the year 
when no sufficiently bright star pairs can be ob- 
served. Importantly, almost all the center-lines 
through the uncertainty ellipses pass very close to 
the pole in the plot. 

The spread between the ACS-derived spin axis 
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Fig. 10. — Histograms of differences between the simulated and observed azimuth (left-hand panel) and ele- 
vation (right-hand panel) of star SA0246574 in Orbit 73. The histograms illustrate the typical determination 
accuracy of star positions from Star Sensor data. 
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Fig. 11. — Histograms of the differences between the simulated and observed orbit- averaged azimuths (left- 
hand panel) and elevations (right-hand panel) of all stars observed by the Star Sensor. 



and that derived by the Star Sensor is between 
0.24°and 0.28°, which roughly agrees with the pre- 
launch measurements of the Star Sensor boresight 
relative to that of IBEX-Lo. From this, we con- 
clude that there is no systematic statistically sig- 
nificant deviation of the IBEX-Lo boresight from 
the direction perpendicular to the ACS-derived 
spin axis. The spread obtained from the analy- 
sis of the Star Sensor data is due to the observing 
conditions and follows a yearly pattern caused by 
the positions of individual stars. This is illustrated 
in Fig. [13] where differences in right ascension and 
declination of the spin axis between the ACS and 
Star Sensor determinations are shown as yearly 
time series based on orbit number. The accuracy 
of the spin axis pointing determined by the Star 
Sensor shows that it could be used for spacecraft 
attitude determination for almost all orbits in case 
of a Star Tracker failure, i.e., the Star Sensor in- 
deed provides good redundancy to the spacecraft 
system. 



6. Summary and conclusions 

We studied the accuracy of the Star Sen- 
sor pointing relative to the IBEX ACS to es- 
tablish the absolute pointing of IBEX-Lo in 
astronomical coordinates during observations 
of neutral interstellar g as flow as p r esente d 



by ^ob iuset al.1 (l2012h: iBzowski et al.l (l2012t): 
fBochsle r et al.l (|2012h : ISaul et al.l (l2012[ ): lLee et al 
(20121 ). We analyzed IBEX ACS and Star Sensor 
data obtained during all orbits over the first two 
years of the mission, for which stars could be 
clearly resolved from background, including orbits 
when the neutral gas observations were carried 
out. 

There are four independent sources of uncer- 
tainty in the determination of the boresight of the 
IBEX-Lo sensor: 

1. Mechanical fabrication and mounting toler- 
ances in the Star Sensor and IBEX-Lo lead 
to a 0.094° azimuth (spin angle) uncertainty 
and a 0.102° elevation angle uncertainty in 
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Fig. 12. — Polar plots of the spin axis pointings in the spacecraft reference system determined from the 
observations of the Star Sensor. The left-hand panel presents separate determinations from all suitable 
spin blocks, the right-hand panel the positions averaged by orbits. The azimuth angles are marked at the 
circumference, the co-latitudes (polar distances) at the co-latitude rings. The z-axis (i.e. the center point) 
corresponds to the spin axis pointing determined by ISOC based on the Stars Tracker/ ACS data. 
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Fig. 13. — Seasonal variation of the differences between the Star Tracker/ ACS and Star Sensor determinations 
of the orbit-averaged spin axis pointings in the equatorial reference system: rectascension (left-hand panel) 
and declination (right-hand panel). Blue lines mark the positions obtained for the orbits from 2009, red lines 
mark the positions obtained for the corresponding orbits from 2010. 



the co-alignment of the Star Sensor and 
IBEX-Lo boresight. 

2. The spin axis is determined throughout each 
orbit from analysis of ACS attitude informa- 
tion is found to be stable to within 0.02° in 
both directions. 

3. Observations of 69 stars and their positions 
were made using the Star Sensor (Fig. [TT]) . 
while the spacecraft ACS was used to de- 
termine the spacecraft spin axis. The mean 
deviation between simulated and observed 
stars in phase angle is equal to —0.0155° ± 
0.017°, while the mean deviation in elevation 



is -0.0179° ± 0.037°. These observations 
show that, on average, there is no statisti- 
cally significant deviation of the Star Sensor 
boresight from the direction perpendicular 
to the ACS-determined spin axis to within 
0.037°, nor in spin angle relative to the spin 
pulse to within 0.017°. 

4. Both the Star Sensor and ACS measure- 
ments are affected by a 1 ms spin period 
granularity, which leads to up to a 0.025° un- 
certainty in spin-phase. 

Table [T] lists each of these measurement uncer- 
tainties and combines them to yield spin angle 
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Table 1: Measurements uncertainties in the IBEX-Lo pointing based on Star Sensor and ACS 





Uncertainty Type 


Spin Angle 
Uncertainty (deg) 


Elevation Angle 
Uncertainty (deg) 


1) 


Star Sensor Mounting Knowledge 


0.094 


0.102 


2) 


Spin axis stability 


0.02 


0.02 


3) 


ACS-determined spin- axis 


0.017 


0.037 


4) 


Spin Period granularity 


<0.025 






Uncertainty in IBEX-Lo Attitude 


0.099 


0.104 




using Star Sensor 








Uncertainty in IBEX-Lo Attitude 


0.101 


0.110 




using ACS 







and elevation angle uncertainties in the IBEX-Lo 
pointing based on instantaneous Star Sensor and 
ACS attitude determination. Uncertainty in co- 
alignment of the Star Sensor and IBEX-Lo bore- 
sights dominates attitude knowledge using both 
the Star Sensor and the ACS. We conclude that no 
systematic correction for boresight geometry needs 
to be introduced in the analysis of neutral inter- 
stellar gas flow observed by IBEX-Lo and the net 
uncertainty in the instantaneous IBEX-Lo point- 
ing is ~ 0.1° in both spin angle and elevation. 

Our analysis shows that Star Sensor informa- 
tion can be used for spacecraft attitude corrections 
during periods when the spacecraft Star Tracker 
(used by the ACS) is disabled by bright objects 
in its field-of-view. In addition, we have demon- 
strated that Star Sensor measurements provide 
independent determination of the spacecraft spin 
axis during most of the orbits with slightly larger 
uncertainties than the ACS-determined spin axis. 
Thus, the Star Sensor can be used for spacecraft 
attitude determination and provides redundancy 
for the Star Tracker in case of any problems. 
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Fig. 14. — Illustration of the sky strip scanned by the Star Sensor during the IBEX rotation. The strip 
is marked in blue; the blue Z" axis marks the rotation axis of the spacecraft, which is also the Z-axis of 
the spacecraft reference system. The black-red Z axis marks the pole of the equatorial reference system, in 
which the stars position are listed in the catalog. The rotation axis is within 1° of the ecliptic plane. 



A. (Appendix) Star Sensor Simulations Program 

The Star Sensor is able to register bright stars, the naked-eye outer planets (Mars, Jupiter and Saturn) 
and the Moon. All of them can be used to determine the pointing of the spin axis and all of them can 
be simulated by the Star Sensor Simulations Program, which is needed to identify the objects observed by 
the Star Sensor with the objects in the reference star catalog and to provide insight into the sky portion 
observed. 

The Star Sensor Simulation Program must identify a list of objects that will enter the FoV for a given 
orientation of the IBEX spin axis at a given time. It is also desirable to determine the objects located in the 
adjacent regions of the sky in the case the true rotation axis differs from the assumed one and the observed 
signal includes the stars that basically should not be visible. Determination of positions of these objects 
on the sky and verification if they are within the visible sky strip at the gi ven moment is done with the 
use of specially developed procedures that use the SPICE software platform ( ActonI [l996[ ) from the NASA 
Navigation and Ancillary Information Facility. 

To determine which stars will be visible for the Star Sensor in a given orbit for a given spin axis pointing 
one must convert the positions of the bright stars from a reference star catalog to the spacecraft frame. From 
the many star catalog systems available for attitude determination, we chose the Smithsonian Astrophysical 
Observatory Catalogue (SAO FK5). The SAO catalog contains almost 260 000 stars down to about 10.0 
visual magnitude. The magnitudes are available only as photographic and/or photo- visual values with an 
accuracy of ~ 0.5 magnitude. The uncertainty of positions is ~ 0.5 arc-sec (Wertz 1978). Because the model 
of the sky background used in the Simulation Program covers the stars up to the magnitude = 6.5, we 
actually used a subset of the SAO catalog that includes only the stars brighter than rriy = 6.5. This was 
done for practical reasons: Star Sensor is not able to resolve weaker stars (actually, the sensitivity limit for 
meaningful applications is about m = 3.5) and Hmiting the source catalog to fewer entries speeds up the 
process of matching the catalog stars with the observed ones. 
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The Star Sensor is able to observe objects located in a strip of the sky between +3.5° and —5.0° from the 
x-y plane of the IBEX spacecraft, which is perpendicular to the rotation axis (see Fig. [T4|) . The rotation 
axis is maintained within a few degrees from the Sun. Thus the field of view of the Sensor changes from 
orbit to orbit. Hence for each simulation for a given date and given pointing of the spin axis the sub-catalog 
is updated for the date of observations and converted from the equatorial coordinate system to the S/C 
coordinate system. The new catalog makes basis for the Star Sensor output signal simulation and scanned 
stars identification. 

Since the Star Sensor is able to resolve the three outer naked-eye planets, the catalog of potentially visible 
stars that is constructed from the subset of the SAO catalog is supplemented with the calculated positions 
of these three planets valid for the time of the simulation. Since the planets are practically point objects and 
their proper motion on sky during the ~ 16 minutes of the duration of one block of Star Sensor observations is 
practically 0, they are treated as stars and included into the Potentially Visible Objects Catalog constructed 
for each simulation. All point-like objects are always observed against a diffuse sky background and since 
the Star Sensor is a pin-hole camera, not an imaging device, the signal registered is a convolved sum of 
contributions from the point sources and the background. The components of the background are the diffuse 
Galactic and inter-galactic light, faint unresolved stars, and the zodiacal light. The first two components of 
the background are practically invariable in time, but the zodiacal light intensity is a function of the angular 
distance of the line of sight from the Sun and thus varies during one orbit. An extensive review of the sky 
background in the optical domain and its components is provided by iLeinert et al.l (|l998[ ) . 

The zodiacal light comes up mostly because of the scattering of sunlight off the interplanetary dust grains. 
Interplanetary dust is distributed symmetrically around the ecliptic plane, so the intensity as seen from the 
Earth orbit is also symmetric around the ecliptic, and because of the optics of the scattering, the intensity 
is a function of the absolute value of difference between the ecliptic longitude of the line of sight and the 
Sun. In the numerical model impl emented for the Star Sensor Simulations Program we adopted the model of 
intensity of the zodiacal light from lLeinert et ah (1998), Table 16, with a correction defined in their Eq. (23). 
The intensity of the zodiacal light is shown in Fig. [151 The FoV of the Star Sensor occupies a 360° band 
perpendicular to the ecliptic, but since the spin axis of the spacecraft is always within a few degrees off the 
Sun, the strip of the sky visible to the Star Sensor covers a very limited band of ecliptic longitudes that differ 
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Fig. 16. — Intensity in the SlOsun units of the diffuse Galactic and starhght as adopted in the model of 
the sky background used in the Star Sensor Simulation Program, shown as function of rectascension and 
declination. 



from the longitude of the Sun by about 90°. Hence the model of the zodiacal light adopted in the Star Sensor 
Simulation Program could be based on a limited range of longitudes differences 75° < |A — AsunI < 105°. 

A nother compon ent of the model of the sky background is the diffuse Galactic light. It was again adopted 
after Leinert et al. 1 ([l998), based on measurements of the sky background from the Pioneer spacecraft per- 
formed beyond the reg:ion occupied by most of the interplanetary dust. We used the data from Table 37 in 
Lein ert et al. (1998'), which are in the SlOsun units. Since the spectrum of the zodiacal light is very similar to 
a slightly reddened solar spectrum, the models of the zodiacal light and of the starlight were compatible and 
since the spectral range of the Star Sensor PMT corresponds to the visible spectrum, it reflected relatively 
well the sensitivity. The data were on a 10° by 10° grid in the equatorial coordinates and covered stars up to 
the visual magnitude m = 6.5. Hence all the brighter stars, even though they could not be clearly resolved 
by the Star Sensor, had to be included in the background model directly from the star catalog. A plot of 
the background from the Galactic diffuse light and unresolved starlight is shown in Fig. [161 

The simulated signal of the Star Sensor from a point source is formed as illustrated in Fig. [T3 Shown is 
the Star Sensor aperture scanning the Visible Strip of the sky. The signal from a point source in the absence 
of the background is shown below the aperture as a double pulse. The output pulses can be approximated as 
triangles with the FWHM determined by geometrical dimensions of the collimator. The relation between the 
star magnitude m and the output voltage U (shown in the figure) is logarithmic. The Star Sensor Simulations 
Program simulates also the random component in the flux, calculating it as a Gaussian noise component 
of the voltage. The relations between the positions of the peaks and their angular separation and the 
coordinates of the scanned star in the S/C reference system are repeated in the figure after Eqs. ([Il[3|). The 
signal observed is a superposition of the signals from all sources. In order to simulate the diffuse background 
an d the stars, w e divid ed the area of the Visible Strip into ~ 800 000 pixels, adopting a tessellation suggested 
bv iGorski et al.l (2005) in their Hierarchical Equal Area Latitude Pixelization (HEALPix) system, mapped 
at the spacecraft coordinates. The equal- area pixels are approximately rectangular in shape and arranged 
in latitudinal rings; the angular distances between the centers of the neighboring pixels are equal to each 
other. These features of this sky tessellation system made it well suitable for our purposes. We mapped 
(separately) the models of the zodiacal light and of the diffuse starlight for a given time and spin axis pointing 
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Fig. 17. — Schematic view of scanning of the sky by the Star Sensor. The strip of the sky visible to the 
Star Sensor is marked in blue. Shown is one bright star within this Visible Strip (black dots), being scanned 
by the split-V aperture of the Star Sensor. Below the Visible Strip the simulated signal is shown. The star 
being scanned makes the two peaks of the output voltage. 

at the HEALPix pixels and then calculated the brightness of each pixel. The model signal from the diffuse 
background is a superposition of the signals from virtual point sources, located at the centers of the pixels. 
The signals from the catalog stars and planets (if present in the FoV) are then added to form the output 
model signal. An example of the fidelity of the model signal as compared with the actual data is shown in 
Fig. El 

In anticipation of an insufficient number of well resolved stars against background, we prepared for deter- 
mination of the spin axis pointing from observations of the Moon. The Moon crosses the FoV approximately 
once per fortnight and it is so bright that the Star Sensor has to be switched into a special Moon Mode to 
prevent the signal from being saturated. The voltage on the PMT in the Moon Mode is reduced, but at a 
cost of reducing the signals from stars below the detectability threshold. The Moon is a fast-moving object 
and its angular size at IBEX locations varies from 0.25° to 4°. Since it is always observed close to the 
half- moon phase, the Simulation Program must take into account its characteristic shape and the fact that 
the center of light is offset from the center of mass position on sky. The Simulation Program projects the 
expected shape of the Moon body at the HEALPix pixels and is able to simulate a realistic Moon signal, but 
because of the very complex nature of the signal and its interpretation on one hand and because point objects 
turned out to be sufficient to fulfill the task of the Star Sensor on the other hand, this mode of operation of 
the Star Sensor Simulations Program has never been exercised to interpret the Star Sensor observations. 

B. (APPENDIX) Spin axis pointing in the spacecraft frame for individual orbits 

Determinations of IBEX spin axis from all suitable spin blocks in the telemetry obtained from Star Sensor 
data were performed in the spacecraft polar reference system. A few of them are shown in Fig. [181 The 
Z-axis coincides with the spin axis pointing determined by the ISOC from ACS data. 
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